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We have analyzed a dissipative two level quantum system (TLS) which is continuously and si- 
multaneously irradiated by a high and low frequency excitation. The interaction of the TLS with a 
high frequency excitation is considered in the frame of the dressed state approach. A linear response 
of the coupled TLS and corresponding photon field system to a signal whose frequency is of the 
order of the Rabi frequency is found. The response exhibits undamped low frequency oscillations, 
whose amplitude has a clear resonance at the Rabi frequency with the width being dependent on 
the damping rates of the TLS. The method can be useful for low-frequency Rabi spectroscopy in 
various physical systems described by a two-level Hamiltonian, such as nuclei spins in NMR, double 
well quantum dots, superconducting fiux and charge qubits, etc. The application of the method to 
a superconducting flux qubit and to the detection of NMR is considered in detail. 
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I. INTRODUCTION 

It is well known that under resonant irradiation a quan- 
tum two-level system (TLS) can undergo coherent (Rabi) 
oscillations. The frequency of these oscillations is pro- 
portional to the amplitude of the resonant fieldi and is 
much lower than the gap frequency of the TLS. This ef- 
fect is widely used in molecular beam spectroscopy^, and 
in quantum optics^. During the last few years it has 
been prove experimentally that Rabi spectroscopy can 
serve as a valuable tool for the determination of relax- 
ation times in solid-state quantum mechanical two-level 
systems, qubits, which are to be used for quantum infor- 
mation processings^. Normally, these systems are strongly 
coupled to the environment, which results in the fast 
damping of the Rabi oscillations. It prevents the use of 
conventional continuous measurements schemes for their 
detection, although special schemes for the detection of 
coherent oscillations through a weak continuous mea- 
surement of a TLS were proposed in^'*''''. Consequently, 
Rabi oscillations are measured by using a pulse tech- 
nique which harnesses the statistics of switching events 
of the occupation probability between two energy levels 
with the excitation and readout being taken at the gap 
frequency of the TLS, which normally lies in the GHz 



Recent successful development of the method of low 
frequency characterization of flux qubits by weak contin- 
uous measurements in the radio frequency domain (see 
review paper^^ and references therein) allowed the first 
spectroscopic monitoring of Rabi oscillations with a low- 
frequency tank circuit which had been tuned to the Rabi 
frequency of the flux qubit These experiments stimu- 
lated theorists to study different methods for the detec- 
tion of Rabi oscillations involving low frequency (com- 
pared to the energy gap between the two levels) electronic 
circuitryiiOSiiejI. 

One of the methods for the detection of Rabi resonance 
at low frequencies has been suggested in paper a^^'^^ . The 



method involves irradiating a TLS continuously and si- 
multaneously by two external sources. The first source 
with frequency luq, which is close to the energy gap be- 
tween the two levels, excites the low-frequency Rabi os- 
cillations. Normally, Rabi oscillations are damped out 
with a rate, which depends on how strongly the system 
is coupled to the environment. However, if a second low- 
frequency source is applied simultaneously to the TLS it 
responds with persistent low frequency oscillations. The 
amplitude of these low-frequency oscillations has a reso- 
nance at the Rabi frequency and its width is dependent 
on the damping rates of the system. In paper o^^'^^ we 
analyzed the effect in the frame of Bloch equations for the 
quantities (o'z), (cy), (cx), where the brackets denote 
the averaging over the environmental degrees of freedom. 
Two external sources at high and low frequency were in- 
corporated in the structure of Bloch equations from the 
very beginning. We showed analytically as well as by 
direct computer simulations of the Bloch equations that 
the quantities (cr^), (cry), (cx) exhibit undamped oscil- 
lations with a resonance at the Rabi frequency. 

The present paper differs fro m^^i^^ in that here we 
study the problem within a dressed state approach, which 
is well known in quantum optics^''. This approach has a 
great advantage over the method used ini^. In that paper 
the effect appeared purely mathematically as the solution 
of the nonlinear Bloch equations for the averaged quan- 
tities, which are not quantum in their nature. It was not 
clear whether the effect had a quantum nature or it was 
a consequence of nonlinear structure of Bloch equations. 

By using the dressed state approach we show here that 
the effect has definitely a quantum nature. From the 
Heisenberg quantum equations of motion we derive the 
Bloch like equations for the elements of the reduced den- 
sity matrix and find the low frequency susceptibilities of 
the TLS coupled to a photon field. We show that, as 
distinct from the cases in quantum optics, in dissipative 
solid state TLS there exists an interaction which can in- 
duce transitions between the dressed Rabi levels. These 
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transitions result in the low frequency response of the sys- 
tem, with the resonance being at the Rabi frequency. In 
addition, the dressed state approach allowed us to find 
important contribution to slow dynamics which comes 
from the transitions between the levels which are spaced 
apart by the gap frequency luq. This contribution is de- 
scribed by the low frequency dynamics of k matrix ele- 
ments (see below), and it cannot be obtained within the 
approach used ini^. Another advantage of the dressed 
state approach is that the expressions we obtain for low 
frequency susceptibility for two level system have a rather 
general nature and can be applied for the investigation 
of specific physical systems. As an illustration of our ap- 
proach in the paper we study in detail two rather different 
physical systems: the flux qubit and real spin under NMR 
conditions. Our approach allows us also to obtain some 
interesting results which could not be obtained within 
the method used ini^. For example, we show that the 
population of Rabi levels becomes inverted under high 
frequency irradiation, that can lead to the heating of the 
low frequency source as was shown in Section 7A for the 
resonant tank of the flux qubit. 



II. INTERACTION OF A TLS WITH A LASER 
FIELD IN THE FRAME OF DRESSED STATES. 
A BRIEF OVERVIEW 

We start with the Hamiltonian of a TLS subjected to 
a high frequency field: 

H ^ ja^ + ^a, + hLUo{a+ a + 1/2) + H,„t (1) 

Here the first two terms describe an isolated TLS, 
which can be used to model a great variety of situa- 
tions in physics and chemistry: from a spin-(l/2) particle 
in a magnetic field to superconducting flux and charge 
qubits^i^. In order to be exact we use the TLS in (jT]) 
to describe a double-well system where only the ground 
states of the two wells are occupied, with A being the 
energy splitting of a symmetric (e = 0) TLS due to quan- 
tum tunnelling between two wells. The quantity e is the 
bias, the external energy parameter which makes the sys- 
tem asymmetric. The third term in ^ is the Hamilto- 
nian of the laser mode, a"*" and a being the creation and 
annihilation operators. The last term in ([T]) describes the 
interaction of the TLS with a laser field. This interaction 
modulates the energy asymmetry between the two wells: 



The paper is organized as follows. In Section II, for 
the sake of the completeness we give a brief overview 
of a TLS interacting with a one mode laser field which 
is tuned to the gap of the TLS. The structure of the 
energy levels of the global system (TLS and associated 
laser field) consists of manifolds where the spacing be- 
tween the two levels in a given manifold is equal to the 
Rabi frequency. We write down the wave functions for 
these two levels and calculate the transition amplitudes 
between them which result from the low frequency ex- 
citation. In Section III we define the density matrix in 
an uncoupled basis and write down the phenomenolog- 
ical rate equations for the elements of the density ma- 
trix. In Section IV the elements of the reduced density 
matrix (the density matrix traced over the photon num- 
ber N) are defined, and the Bloch like rate equations 
for the elements of the reduced density matrix are de- 
rived. For the case of small high frequency detuning the 
steady state solutions to these equations are found. It is 
shown that under high frequency irradiation the popula- 
tion of the Rabi levels becomes inverted. In Section V 
the Bloch equations for the reduced density matrix are 
modified to include the low frequency excitation. The 
linear low frequency susceptibilities for the response of 
the coupled TLS and associated photon field system to 
a low frequency excitation are found in Section VI both 
for arbitrary and small high frequency detuning. Two 
important applications of the method (superconducting 
flux qubit and NMR spins) are considered in detail in 
Section VII. 



H,nt-- -^a,F{a+ + a) (2) 

The Hamiltonian of the TLS in ^ is written in the 
localized state basis, i.e. in the basis of states localized 
in each well. In terms of the eigenstates basis, which we 
denote by upper-case subscripts for the Pauli matrices 
<7Xi o'Y, and crz, Hamiltonian ([1]) reads as: 

H = ^az + tiiooia+a+ 1/2) + H,„t (3) 

where = y/ A^ -I- is the gap between the two energy 
states and 

First, we consider a noninteracting system which involves 
the TLS and a laser field. The system is described by the 
Hamiltonian: 

Ho ^^az + hLUoia+a + 1/2) (5) 

We denote the ground state and excited state wave func- 
tions of the TLS as \a) and \b), respectively, where they 
have properties: crz\a) — — |a), az\b) = |6), (Jx\a) = 
cxlb) = I a). The eigenfunctions of the photon field are 
|iV): a+\N) = ^/WTT\N + 1), and a\N) = VN\N - 
1). We denote the eigenfunctions of the noninteracting 
TLS and associated photon system as a tensor product 
\a,N) = \a) ® |iV), and \b,N) = \b) (g) \N). Up to a 
constant term the energies of these states are: 

Ho\a,N)=(-^ + nLUoN\\a,N) (6) 
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I a, N+2> 
|i',N+l> 

I a, N+l> 
lb,N> 



|a,N> 
|fo,N-l> 



In quantum optics there is no "longitudinal" interaction 
between an atomic spin and photon field which is propor- 
tional to az- It is the presence of this bias interaction 
in a dissipative TLS that leads to some effects which are 
not observed in quantum optics. 

It should be considered now how these levels are mod- 
ified due to the interaction with ([4]). For a pair of these 
closed spacing levels within a given manifold, \a,N), 
\b, N~l) interaction (|4]) causes a transition between them 
with the amplitude 



AF 



/ \ r , 

{a,N\H.M\b,N-l)^—VN 



(8) 



Therefore, within 2D Hilbert space the wave functions 
|a, N) and \b, N — 1) are mixed to give new wave functions 
for the dressed states \l,N), and \2,N): 



smf a 



N) +cos9\b,N -1) 



(9) 



|a,N-l> 
I b, N-2> 

I a, N-2> 
I b, N-3> 



\2,N) = -cos6'|a,iV) +sin6'|6,iV- 1) (10) 

where we define the state with higher energy as |1, TV). 

The form of Eqs. ^ and ifTO]) ensures the normaliza- 
tion and orthogonality of the wave functions |1,-/V) and 
\2,N), which are eigenfunctions of Hamiltonian ([3]). Ac- 
cordingly, the uncoupled states \a,N), \b,N — 1) can be 
expressed in terms of the dressed states |1, A^), |2, A^): 



\a,N) = sin 6*11, A^) -cos6l|2,A^) 



(11) 



FIG. 1: (Color online). Levels of the noninteracting TLS and 
associated photon field system. The level spacing in a given 
manifold is hS. The spacing between neighboring manifolds 
is huiQ- The red (blue) arrows show the transitions induced 
by Hint from the state \b,N - 1) (|a, A>). 



Ho\b,N)^i^^ + nwoNj\b,N) (7) 

For a photon frequency luq close to the TLS frequency 
A^/S, and a small detuning S = ujq — A^/h ^ ujq, Ag/ft, 
where for definitiveness we assume i5 > 0, it is seen from 
^ and ([7]) that the energies of the states \a, N + 1) and 
\b, N) are close to each other: Ea,N+i — E^^m = fiS. The 
same is true for the pairs of states |a. A') and \b, N — 1), 
\a, N + 2) and \b, N + 1), and so on. Therefore, the en- 
ergy levels of the system under concideration is a ladder 
of pairs of manifolds which are specified by the photon 
number A^ (see Fig. [T]). Every manifold is parameterized 
by a pair of states with a small spacing between them, 
hS, and the distance between neighboring manifolds is 
equal to the photon energy, Hluq. 

This ladder of manifolds is quite similar to the one 
for atom-field interactions^*'. However, a principal differ- 
ence is the structure of the interaction Hamiltonian Q. 



|6,Ar- 1) = cos6l|l, A^) +sin6l|2, A^) (12) 

By using standard quantum mechanical techniques we 
can find the eigenenergies and the angle 9: 

E± = \ {E\a,N) + E\b,N-i)) ± \hnR (13) 

where the upper (lower) sign corresponds to |1,A^) 
(1 2, AT)). The quantity in ^ is the Rabi frequency^S 

= ^j5^ + nl (14) 

where Vli — AF/hAg, and ^/N is incorporated in the 
high frequency amplitude F. 

For the angle 9 we obtain tan26' = —fli/S, where 
< 20 < vr, so that cos 20 = -5/VLr, cos9 = 

72(l-l4)'^''^-dsin0=-l5(l + ^)'^'. 

When the interaction is switched off {F — > 0) then, 
as would be expected, the state |1,A^) tends to \a,N), 
and the state |2, A^) tends to \b,N — 1). Therefore, by 
taking into account the interaction between the TLS and 
the photon field the level structure of a given manifold 
looks like that shown in Fig. [21 The interaction results in 
an increase of the energy gap between the states \a^N) 
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1,N) 



|a,N> 



,N-1) 



|2,N) 



FIG. 2: The level structure of a dressed manifold. The un- 
coupled states \a,N) and \b,N ~ 1) are transformed by the 
interaction between the TLS and photon field into the two 
dressed states lljA^) and |2,7V). 



III. RATE EQUATIONS FOR THE DENSITY 
MATRIX IN UNCOUPLED BASIS 



The elements of the density matrix a in the basis of the 
uncoupled levels |a, iV) and \b,N) arc defined as follows: 



NN' 
aa 

NN' 
'bb 

NN' 
'ab 



{a,N\a\a,N') 
{b,N\a\b,N') 
{a,N\a\b,N') 



(17) 



Since the spontaneous transitions between the levels \a) 
and 1 6), which lead to their decay do not change the pho- 
ton numbers the rate equations for the density matrix a 
can be written in the following form: 



and \b, N — 1). Consequently, it is said that these states 
are dressed by the interaction. From now on we will call 
these two nearby dressed states Rabi levels. 

Up to now the picture has been quite similar to that 
known from atom-photon interactions^. However, a 
drastic difference appears if we consider the excitation 
of the dressed levels |1, A^) and |2, iV) by a signal the fre- 
quency of which is comparable with Rabi frequency Qr. 
Such a low frequency signal cannot change the number 
N of high frequency photons, therefore, in quantum op- 
tics the transition between these two states is not allowed 
since, in the language of quantum optics, the atom dipole 
operator is transversal: it is proportional to ax- There- 
fore, the atom dipole operator connects only the levels, 
say |a, N) and \b, N), which belong to the different man- 
ifolds. 

In the TLS the transitions between nearby states \ l;N) 
and 1 2, TV) are caused by the "longitudinal" term in 
the interaction Hamiltonian, which is proportional to 
az- Let us assume that system ([T]) additionally inter- 
acts with the low frequency signal G cos cut, so that the 
term — a zG cos cut is added to Hamiltonian ([1]), where the 
frequency lu is of the same order of magnitude as the 
Rabi frequency Q^. In the eigenbasis this low frequency 
Hamiltonian is transformed as follows: 



"'"bb 
dt 



da. 



NN' 

aa 



dt 



-^{^bN - ^bN')(^ab ~^bb ii+C^aa T 

(18) 



* /77 p \„NN' „NN'-r I „NN'-r 

--^[j^aN - -t^aN'jCr^f, - a ll+Cr^f, I4 



dt 



{EaN — EbN')cra 



NN' 
b 



„NN't^ 

a„h i L, 



(19) 



(20) 



where is the transition rate from the state |&, iV) to 
state |a, iV) (relaxation rate), F| is the transition rate 
from the state |a, TV) to state \b,N) (excitation rate), 
and the quantity F,^ is the rate of decoherence. For equi- 
librium conditions the relaxation and excitation rates are 
related by the detailed balance law: 



Tt = Ti exp - 



(21) 



From (PT|) we obtain 



LF 



A, 



■ax 



az I G cos ojt 



(15) 



F_Ti 



tanh 



A, 



2fcsT 



(22) 



The transitions between the Rabi levels |1, TV) and |2, TV) 
are caused by the second term in the brackets of P^ : 



(l,TV|i/fJ|2,TV) 



eGVLi 



R 



(16) 



In what follows we show that these transitions result in 
the undamped low frequency oscillations of the popula- 
tions of initial levels \a) and \b) which can be detected by 
using appropriate electronic circuitry. 



where we define F_ = F| — F| and the longitudinal re- 
laxation time T]~^ = F| -t- F|. 

The equations (|18p. and (|^D|) can be written in 

the operator form: 



(23) 



where the Hamiltonian of the uncoupled system ffp is 
given in The operator L is defined by its matrix 
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elements which follows from ^8^, (|T9]) and pO)) : 

L^-r^ \bNi){bNi\a\bN2){bN2\ 

Ni,N2 

+ r| ^ |5A^i)(a7Vi|(7|a7V2)(6iV2| 

-Ti ^ |aiVi)(aiVi|(7|aiV2)(a7V2| 
N1.N2 

+ ri \aNi){bNi\a\bN2){aN2\ 

Ni,N2 

-r^ ^ \aNi){aNi\a\bN2){bN2\ 

Ni,N2 

-r^ ^ \bNi){bNi\a\aN2){aN2\ (24) 

Ni,N2 

Below we assume that the interaction between the TLS 
and laser field does not influence the F's rates and the 
corresponding rate equations. As was shown i n^3,24 ^ ^j^jg 
assumption is valid for relative weak driving, sufficient 
short correlation time of the heat bath Tc, and in large 
temperature hmit: F « A^, S/tc, ksT—. 



IV. BLOCK TYPE EQUATIONS FOR THE 
REDUCED DENSITY MATRIX 



Equation ((23)) can be generalized to include the inter- 
action between the TLS and laser field: 



da 
It 



-[H,a] +L 



(25) 



where the Hamiltonian H is given in ^ with Hint in Q • 
Following on, we define the reduced density matrix for 
the two level coupled system by tracing over the photon 
number N: 



N 
N 

pi2 = 5](l,7V|a|2,iV) 

N 

P21 -5I<2,7V|a|l,iV) 



(26) 



AT 



^n = ^(l,iV|a|l,7V-l) 

N 

«;22=^(2,iV|a|2,7V-l) 

N 

^i2=^(l,iV|a|2,7V-l) 

N 

«;2i=5](2,iV|a|l,7V-l) 



(27) 



It is convenient to write the rate equations in terms of 
the new variables: p = pn - P22, P+ = P12 + P21, P- = 
P12 - P21, = Kii + K22, K = Kii- K22, n+ = K12 + H21, 
K_ = K12 — ^21- It is not difficult to show that the total 



population p+ — pn + P22 is constant: 



dt 



with the 



normalization condition = 1. 

By taking into account the matrix elements of L in the 
dressed state basis (see Appendix A) , we obtain the rate 
equations for the elements of the reduced density matrix: 



dp . 



Bp, 



r_ cos 29 



dp 
~dt 



= -iflnp- + Bp - A2P+ + r_ sin 29 



dp- 



dK+ 



-VjJqK 



dn, 



— — —iujQK — AiK + Bk+ + K^r_ cos 29 



(28) 

(29) 

(30) 
(31) 

(32) 



— ^ — zcjoK-i-— if2flK_-|-i3K— yi2K+-|-K r_ sin 2w (33) 



where 



dK- 
~dt 



A2 



B = 



— cos^ 29 + sin^ 29 



sin^ 261 + cos^ 261 



sin 26* cos 261 



(34) 

(35) 
(36) 
(37) 



JV 



The p's elements of the density matrix describe the 
transitions between the Rabi levels, which are usually 
accounted for by a so called rotating wave approxima- 
tion (RWA). The /c's elements describe the transitions 
between neighboring manifolds spaced apart by the en- 
ergy hujQ. In the above equations it can be seen that these 
two types of transitions are completely independent: the 
equations for the p matrix are uncoupled from those for 
the K matrix. 

If the damping is absent (aU T's in (gSl), ([M]), (001) are 
equal to zero) the quantity p is constant, and p+ and p_ 
oscillate with the Rabi frequency Qr. However, in the 
presence of damping these oscillations rapidly decay to 
their steady state values. Here it is instructive to consider 
two limiting cases. 
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A. Bloch equations in the absence of a high 
frequency excitation 

Let the power of high frequency photon field tends 
to zero. In this case fti ^ 0, ojq 0, and therefore, 
sin26' — + 0, cos 20 — > ~1, and A^/h. From Eqs. 

I (EH), dSOl) we get: 



dp 1 ^ 



dp+ 
dt 



(38) 



(39) 



C. Steady state solution for the density matrix 



The steady state solution (^ = ~ = 0) for 



Eqs. dill), jSO]) is as follows 



r2 

Ti 



r_ cos 26* 



(0) 

P+ 



r2 

Ti 



-r_ sin 261 



JO) 



.^^fl (0) 



(44) 



(45) 



(46) 



~dr 



(40) 



This result is obvious: the ladder of manifolds is reduced 
to just two levels |a) and |6). The diagonal part of the 
density matrix relaxes with a rate of 1/Ti to its steady 
state value p^^^^ — — TiF^, while the offdiagonal elements 
p± exhibit the damped oscillations with a rate and 
frequency Ae/ft. 



It is interesting to note that under high frequency ir- 
radiation the population of the Rabi levels becomes in- 
verted. In (UH), the quantity p^^\ which is the dif- 
ference of the populations between higher and lower 
Rabi levels, becomes positive, since for 5 > we have 
cos 20 = < 0, and always F_ < 0. 

For the case when the high frequency detuning 6 is 
small compared to the Rabi frequency, at zero detuning 
{6 ^ ill), we have sin26' 1 and cos 20 — (J/ili, and 
we get from Eqs. dH]), (gS]), and (|46)) : 



B. Bloch equations for zero high frequency 
detuning 5. 

Consider the case when the high frequency detuning is 
zero ((5 = 0). In this limit sin 20 1, cos 20 — > and we 
get from the Eqs. 1^, (EH), (P| : 



dp 
'di 



dp+ 1 

-df ^ -'""'P- t/^ 



dp- 
~dt 



= —iftip+ — Fyp_ 



(41) 



(42) 



(43) 



It is seen that the high frequency excitation drastically 
changes the behavior of the density matrix. Here the pop- 
ulation p decays with the decoherence rate F^. This is 
due to the fact that the population of the level, say, \1,N) 
can only be changed as a result of spontaneous transitions 
to the levels |1, TV — 1) and |2, iV — 1) of the neighboring 
manifold. The offdiagonal quantities p+ and p- undergo 

the damped oscillations with a rate i P™" 
vided that > ^ - F^ 



Ti 



.(0) 



F-Fy 

"11 Ti 



F_f2i 
i 



(47) 



(48) 



(49) 



As seen from (|47p p^ ' — + as (5 tends to zero. This 
causes the equalization of the population of the two levels 
(pii = P22 = 5) when the high frequency is in exact 
resonance with the energy gap of the TLS. 

The steady state solutions of Eqs. 
and are equal to zero: = 0, k'^^^ 



0, = 0, 



,(0) 

,(0) 
-^21 



■ 0. This implies nf-^^ — 0, 
0. 



,(0) 



— 0, K^2^ — 0, and 



V. EXCITATION OF RABI LEVELS BY A LOW 
FREQUENCY SIGNAL 

In this section we find the response of the TLS coupled 
to a photon field to an external signal the frequency of 
which is of the order of the Rabi frequency Qr. The 
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operator equation for the density matrix a is similar to 
13: 



da 



(50) 



where the Hamiltonian Hj^^ is given in p^ . 

Since the low frequency signal cannot change the pho- 
ton number N, the transitions between Rabi levels |1, iV) 
and |2, iV) can be induced only by the second term in the 
low frequency Hamiltonian (|15|) . The equations for the 
p's and the k's are obtained in the same way as used 
for Eqs. pS)) - The only difference is the appear- 

ance of low frequency terms in the right hand side of 
these equations. Therefore, by taking into account the 
low frequency excitation, we get the following Bloch like 
equations for the reduced density matrix: 



dn- 



— ig-z cos Lut [p+ sin 29 + p cos 29] 

+ igi cos ojt [k+ cos 29 ~ n sin 29] (57) 
The equation for k+ remains unchanged (Eq. (j3ip ). 



Therefore, the low frequency excitation couples the p and 
K elements of the reduced density matrix. Due to this 
coupling the k, elements acquire the low frequency part 
which thus should be included in the standard scheme of 
RWA. 



VI. LOW FREQUENCY LINEAR 
SUSCEPTIBILITIES FOR THE TLS 



dp 
'dt 



-Aip + -Bp+ — p- (igi sin 29 cosujt) 



ig2 cos Ljt cos^ 9 (ki2 — 



r_cos26' (51) 



—j^ = —ig2C0Sujtsm 6* (k2i — kJ^) (52) 



= —iQfip^ + Bp — A2P+ + p- {igi cos 261 cos wi) 

+ ig2C0sujt[{K22 — K22) + {i^ti ^ '^ii) cos 20+ 
(k2i + Kjt, - «2i - ^12) sin 29] + r_ sin 29 (53) 



= —ifljiP^—T^pp^+igi (p+ cos 29 — psin 29) cos ivt 

— ig2 cos [sin 20 {ki2 + '*i2 + '*2i + ^21) 

- (k11 + - {k22 + Kt2) COS 29] (54) 

where gi = 2eG/hA^, and 52 = 2AG/hA^. 



It is evident that the above equations exhibit oscilla- 
tory solutions in the presence of damping. For a small 
amplitude low frequency signal, G the time dependent 
solution for Eqs. (HI]), and ^ can be ob- 

tained by using a perturbation method where the small 
time-dependent corrections to the steady state values 
are: p{t) ^ p(°) + p(i)(t), p+{t) = pf + p^l\t) , and 

P-{t) — p^'' + p^}_\t). By doing this we may neglect all 
k's on the right hand side of these equations, since the 
steady state values for the k's are zero. Therefore, in 
this approximation the equations for the time dependent 
corrections to the p's are decoupled from those to the k's: 



dpW 
dt 



~Aip 



(1) 



BP^^ 



dp 



(1) 



dt 



-inRp^l'>+Bp^^^- 



p^°\ig sin29 cos ojt) (58) 



-A2P4^^-l-p^'' {ig COS 29 COS Lot) 



(59) 



— = —iujQfi — AiK + Bkj^ + K^T^ cos 29 
dt 

-\- ig2 cosujt [p_|_ — p_ cos 20] — igiK^ coscji sin20 (55) 



ifljiK^ + Bk — ^2^?+ + K^r_ sin 29 



dK+ 
dt 

— ig2 cos Lot [p_ sin 20 + p] + igin^ cos cot cos29 (56) 



dpW 
dt 



ig (p^"^ cos 20 - p'-°^ sin 20^ cos ut (60) 



where p(°) , p^-* , and pf^*'' are the steady state values given 
in (m, gSD, and 1^. 

From these equations it is not difficult to find the linear 
susceptibilities of the system {xp{^) = p{^)IG{uo), etc.): 



J 



Xp(w) 



2eVL, 



sin 20 



{i^ + r^) ( + ^ ) + r!| 



02 

-^B cos 29 



(61) 



8 



{iuj + T^) [iuj + — ] - (iuj + Ai) 



XP- i^) - -^ n/^^T^A pT cos 2^? f zc. + 1-) (tuj + 2r^) 



D{Lo)hTlA. 



Ti 



where 



D{uj) = {iuj + T^f [iLO + ^]+{iijj + Ai) nj 



(62) 
(63) 

(64) 



r 



Eqs. ([6T|) . ([62|l . and ((63|) give the response of the ecu- as follows: 
pled system (TLS and associated photon field) to a low 
frequency signal, which excites transitions between the 
Rabi levels. 

From ([5^ . and ([55)1 the linear susceptibilities for 
the case of small high frequency detuning, sin26' 1, 
cos 26* -6/ni, Ai r^, A2 1/Ti, and B 

— jIj- — , can also be obtained. 



dK+ 

IT 



-VjJqK 



(69) 



= —iujoK — AiK + Bk^ + K^r_ COS 29 



i^ficosuot (70) 



2£»i 



2s pT 



where 



d{uj) = (iuj + {iuj + r^) + nj 



(65) 



(66) 



(67) 



(68) 



dn^ + ■ /I 

— ; — = —iluoKj^ — iilfiK^ + Bk — A2H+ + K r_ Sl\i26 
dt 

-i^f2COfiLdt (71) 
— ^ = —iujQK.^ — ifliiKjf- — r^K_ — ^^■^s coswt (72) 



where 



-pL"^cos2^? 



/2= (pL°^sin20 + p(": 



(73) 



(74) 



and p^i"'' is given by (|15|) . The resonance nature of the 
response is evident from ((68)) . 

The equations for the time dependent k's in the first 
order in g are coupled to p's via their steady state values 



J 



/3 = (^p^^'> sin 261 + cos 26*) (75) 



The low frequency linear susceptibilities for the k's (xt^i'-^) — '*('^)/G(w), etc.) are as follows: Xn+i^) = 0, 



Xk- (^) 



XK(wj 



X«+('^) 



A 

A 



(-i/3 [i(c^ + c^o) + Ai] [i{u; + ujo) + A2] + fiBflH + ihB^ ~ /af^fl + uj^) + Ai]) (76) 
(z/i [i(cj + cjo) + v42] [z(c^ + c^o) + r^] - ihB \i{io + cjo) + T^] - hB'^R + i/s^^) (77) 



{-if 2 [i(w + Wo) + ^1] [iiuj + too) + r^] + ifiB [i{uj + too) + T^] - f^nR [i(w + uq) + Ai]) (78) 
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where 



Do{lj) = [i{Lj + cjo) + Tipf + ^o) + jr^ 



(79) 



By assuming that the gap frequency ujq is much more 
than oj and the rates F's, we obtain from (|75|) . ([771) . 
(|78p the corresponding low frequency susceptibihties in 
the second order of the inverse frequency loq: 



-h 



A 



(80) 



-/2 



{Alh + hB) 



(81) 



(82) 



VII. THE APPLICATIONS 

In possible applications of the method that we pro- 
pose here the quantities to be measured are the averages 
of the Pauli spin operators (crx)j (fy), (cz), and their 
time derivatives '^^^^^ , , ^^^^ ■ By using the defini- 
tions of the density matrix ((26)) and (p7|) , and the dressed 
states © and (fTU]) we obtain from the direct application 
of Eqs. (fTTS]) and (fTTB)) (see Appendix B): 



(az) ^ pit) cos 20 + p+{t) sin20 



(83) 



(ax) =sm29Re[K{t)]~cos2dRe[K+{t)]-Re[K^{t)] (84) 



(cry) = -sin26'/m[K(t)] + cos 26'/m[K+(t)] + Im[K^{t)] 

(85) 

For a small amplitude low frequency excitation G cos tot 
we calculate, below, the time derivatives of (|83p . ([84)) . 
and (|55|) with the help of equations (|55|). ([55)) . and 
for the p's and ([701), dZH), and dZH) for the k's: 



= ujo sm29Im[K(t)] - sin26'i?e[K(i)] 



+ cos 2eRe [k+ [t]] - {u;o + ^ b. cos 29) Im [k_ [t)] 
- {ojocos29 + Ur) Im[K+{t)] + V^Re[K^{t)] (86) 



dt 



= LJQ sin26'i?e[K(t)] -f sin26'/m[K(i)] 



- r^cos26'/m[K+(t)] - [ujo + ^rcos29) Re[n^[t)\ 
'[ijjQ cos 20 + riji) Re[K4^{t)]—T ^Im[n^{t)\—g2fz cosuot 



(87) 



- -;i(az)-zl]«sin20p(°^-ir!flSin20pL'^(O+r- 

(88) 

where {uz) in the right hand side of (|88p is written using 
the first approximation only: 

(ffz) = P^"^ cos 261+p^"^ sin 2e'+/9(^) (i) cos 29+p'^l^ (t) sin 26* 

(89) 

The quantities p''^\t), p^^\t), and p^:l\t) can 
be expressed in terms of the real (%') and imagi- 
nary (x") parts of their corresponding low frequency 
susceptibilities: P^^H^) = C (xj, cos wt — x" sin cjt) , 

P^+\t) = Cfx' cosu;t-x" sincoiV and pL'^(i) = 



iG (^Xp_ sinLut + x"p_ cosujtj . 



A. Flux qubit 

Our method can be directly applied to a persistent cur- 
rent qubit (flux qubit), which is a superconducting loop 
interrupted by three Josephson junctions— i^^. For these 
qubits the successful experimental implementation of the 
low frequency readout electronics has been demonstrated 
13,30,31 ^ Though the main subject of the paper is the com- 
bined effect on the TLS of the high and low frequencies, 
it is worth noting that the pure low frequency probing of 
the TLS can also be a valuable tool for characterization 
of the two-level systems. In fact, for the flux qubit this 
method has been realized theoretically and experimen- 
tally 

The basis qubit states have opposing persistent cur- 
rents. The operator of the persistent current in the qubit 
loop is Iq = IqUz- In eigenstate basis the average current. 



is: 



A, 



(e(az)-A(cTx)) 



(90) 



This current can be detected by a high quality reso- 
nant tank circuit inductively coupled to the qubit loop 
13,30,31 ^ggg Y\g. 12]). For the flux qubit the bias e is con- 
trolled by an external dc flux $x: £ = £^j(*i>Js:/*I>o — 1/2), 
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FIG. 3: Flux qubit coupled to a high quality resonant tank 
circuit. 

where Ej = ^oIq/2TT, and 3>o — h/2e is the flux quan- 
tum. The qubit is inductively coupled through the mu- 
tual inductance Af = k{LqLTY/'^ , where k is the dimen- 
sionless coupling parameter, and Lq is the inductance of 
the qubit loop, to a high quality resonant tank circuit 
with inductance Ly, capacitance Ct, and quality factor 



Qt- The tank circuit is biased by a low (compared to 
the gap) frequency current = Iq cos tot. In addition, 
a high frequency signal Fcosbj^t (not shown in Fig. [3|) 
tuned to the gap frequency A^/h is applied to the qubit 
loop. This readout circuit has prove to be successful for 
the investigation of the quantum properties of the flux 
qubit l^i^i^. The interaction between the qubit and the 
tank circuit is described by the term G{t) = Mlqlxit), 
where J, is the current flowing in the qubit loop, and It 
is the current in the tank coil. The voltage across the 
tank, V{t) — Vt cos(a;t + x)j is given by—: 

V + -fTV + u;tV = -Muj^^^^ +uj^LTib, (91) 

where 77^ = ujt/Qt: Etnd ujt — {LtCt)~^^^ is the tank 
resonance frequency, which is tuned to the Rabi fre- 
quency {uiT — ^r)- By using and we find: 



= --r^ [p{t) cos 261 + p+{t) sin 261 + iilRTi sm29p-{t)]-^{ujosm20Im[K{t)]-T^ sin20Re[K{t)]+T^Re[K^{ 
at AgTi 

-I- Fy cos26'i?e[K+(i)] - (t^o + ri_RC0s26') lTn[K^{t)] - (wq cos 26* + fin) Im[K+{t)]} 

(92) 

In order to write down the Fourier component of l|92|) we use the relations (/m[fc(i)])^ = ^ {Xk{^^) — X*k{—^)) G{ijj), 
{Re[k{t)])^ = i (xk(u;) + x%{-^)) G(w), etc... Therefore, 



dt 



IqSG{uj 
AeTi 



[Xp(t^) cos 26* + Xp+ (t^) sin 26* + i^RTi sin 2exp 



IqAGjco) 
2A, 



{-iujQ sin 26* [xk(u;) - x*(-w 



Xk-('^) + XL(-'^) -r^sin26l[xK(w) +Xk(-w)] 4-r<^cos26l Xk+ (^) + xL (-^) 



+ i{ujo + nRCOs20) Xk-('^) - xL(-^) +iiujocos20 + nR) Xk+ (t^) - xL (-^) } 



where the expressions for the susceptibilities are given in (ISTj) . ([5^ . (pS)) . ([77)l . and (|75|) . 

Eq. (j9ip in terms of its Fourier components reads: 



V{lu) {ujj, -uj^ + iuj-iT) = -Mlu^ i j + 



iujujrpLTia 



(93) 



(94) 



By taking into account that G{lu) = M Iqlxii-^), where /t(w) = —iV{Lu)/ujLT, we obtain for the low frequency 
detuning ^, and friction, Tt- 



2 2 k'^^TLql^e 



k^LJTLql^A 

2A, 



i^-^ (nnT, sin 20Xp_ (c^) + cos 20x"(^) + sin 20x"+ ('^)) 

(-^0 sin 20 [X'„(^) - x«(-^)] [xL(^) - X'^J-'^)] - sin 20 [x'^{u;) - x"(-^) 



cos 26* 



X"+ (c^) - Xk'+ (-^)] + K cos 20 + ^r) [x',^ (to) - x'.+ (-cc-)] + (c.0 + nR cos 20) [xL (^) - xL (-'^)] ) 



(95) 
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It 



cos 29 



^ (tIrTi sin 20x;_ (^) - cos 20Xp(^) - sin 20x;+ 



OJ 



2A, 



(too sin 20 [x'^(a>) + x"(-^)] + [xL (^) + xL ("w)] " sin 20 [x',{lo) + X^-^)] 
Xk+ (^) + Xk+ (-^)] - (^0 cos 29 + Qr) [x:^ (u;) + x'L^ {-u)] - (c^o + ^r cos 26) [xL (^) + xL (-'^)] ) 



(96) 



From ([55)) . and ([M)) we obtain the voltage amplitude 
Vt, and the phase, X'- ~ cjcj^Lt/qZ-v/^^ + "^^r^, and 
^5X = 'f/t^rT. 

These expressions have two different parts. The terms 
which are proportional to e are due to transitions between 
the Rabi levels. These terms have a resonance at the Rabi 
frequency and they vanish at the optimal point (at e = 
0). Other, nonresonant, terms come from the transitions 
at the frequency between neighboring manifolds. The 
contribution of these terms does not vanish at the optimal 
point. 

By assuming that the gap frequency is much more 
than u) and the rates F's, we obtain at e = with the help 
of the susceptibilities ([5T|) . and (|82p the expressions 
for low frequency detuning ^ and friction Tt'- 

i.2 T t2, ,2 



Similar effects have been observed These experi- 

ments showed an increase in the quality of the tank res- 
onator by approximately a factor of 2. For the estimation 
of Ft we take the parameters of the flux qubit used in^'^: 
A/ft = 1 GHz, Lq = 24 pH, ujt/^tt = 6.284 MHz, the 
tank quality factor Qt = 1850, jt = ^t/Qt — 2.0 x 10*, 
/, = 600 nA, Lt = 0.2/iH, M = 70 pH, F^ = 8 x lO^c-^, 

and ni - ujT- The result is Ft « 7t - 2 x 10"* . For 

5 O.lcjQ we obtain an increase in the quality factor of 
the tank by approximately 15%. Therefore, for the pa- 
rameters we used expression ((98|) gives a relatively weak 
effect. 

As was shown by Ju. Hauss et al. (2007)^^, an addi- 
tional contribution to the heating of the tank, which can 
explain the observed amplitude, appears in the second 
order terms in G (which is beyond our linear approxima- 
tion) at twice the Rabi frequency lo ~ 20. r. The effect 
appears at nonzero high frequency detuning i5 as a result 
of the "negative temperature" of the Rabi levels men- 
tioned before in connection with Ea. (|l^ . 



It 



k^Lql^ 



20fj 

V 



+ 



TiF3 



where 



~(0) 



(98) 



(99) 



At the point of resonance {S = 0) the response of the flux 
qubit is purely inductive. Ft = 7t- It is interesting to 
note that expression ((98| predicts the " cooling down" of 
the irradiated qubit due to its interaction with the tank 
circuit. Since F_Ti is negative (see (1221) 1 the quantity 
Ft at (5 > is less than 7t, i. e. the quality factor of the 
tank is increased. The energy from the qubit transfers 
to the tank: the qubit is cooled, the tank is heated. The 
effect is more pronounced if {S ^ fli). For this case 
expression (|98)) is reduced to: 



Ft — "fT 



k^LJ] 



q^q 



1 ^ A \ 

u, tanh 

^ \2kBTj 



(100) 



B. Nuclear magnetic resonance 

Another application of the method can be made to nu- 
clear magnetic resonance (NMR) . The basic input scheme 
of any NMR device is similar to that shown in Fig. [31 
where a sample with a substance under study replaces the 
flux qubit loop. The common mode of operation of the 
NMR device is to polarize the sample with a relatively 
high static magnetic field Bq and at the same time apply 
a time dependent magnetic field Bi cos toot perpendicular 
to Bq. The amplitude of the excitation signal is rather 
low {Bi Bq) and it is tuned to the NMR resonance fre- 
quency {luq "fBo), where 7 is the gyromagnetic ratio. 
The tank circuit detects the time derivatives (Faraday 
law of induction) of the transversal magnetizations Mx 
and My which oscillate with the high NMR frequency 
7_Bo. In addition, under the excitation Bicoscuot the 
longitudinal magnetization Mz oscillates with the Rabi 
frequency ^Ir « 7^1, and sidebands luq ± Qr appear 
in the transversal components Mx and My- However, 
due to coupling to the environment, which is described 
by the relaxation Ti and dephasing T2 times in Bloch 
equations, the Rabi oscillations undergo fast decay and 
disappear from the output signal. The application to the 
sample of another low frequency signal B^p cos cut, which 
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is tuned to the Rabi frequency (w « rifi) causes a persis- 
tent low frequency oscillation of the magnetization with 
its resonance at the Rabi frequency. 

The case of NMR differs from the previous example in 
that here we deal with a real spin-1/2 particle. There- 
fore, all components (di) {i = X, Y, Z) of the averaged 
spin operator and their time derivatives are accessible for 
the measurements. As distinct from the flux qubit case, 
where the low and high frequency excitations are cou- 
pled only to az (see Eq. [2]) in NMR we may couple the 
low frequency and high frequency fields separately to any 
component of the spin operator. In order to simplify the 
problem we consider the following Hamiltonian for the 
spin-1/2 particle in an external magnetic field: 

H = — - — <7z + —<TzBLFC0sujt+ —axBi coswo* 

(101) 

where we choose the z axis to be along the polarizing 
field Bq, and the sigh in (|10ip is consistent with our con- 
vention az\a) = —\a), crz\b) = \b). 

The high frequency excitation signal Bi cos ujQt, which 
is tuned to the NMR resonance frequency ^Bq, is applied 
along the x axis, and a low frequency signal B^p coswt, 
which is tuned to the Rabi frequency, lo w Vlji, is applied 
along the z axis. It is just the low frequency component 
in (|10ip that causes the transitions between the Rabi 
levels. 

Since the low frequency component is applied only 
along one axis (as distinct from (|15p ). (72 = 0, and equa- 
tions ([?T1) - ([51)l for the p's are decoupled from equations 
(|55|) . and (|57|) for the k's. This leads to a nonzero 

solution for the k's only in the second order in low fre- 
quency amplitude gi. Therefore, in this case we have 
low frequency responses only for the quantities (az) and 
{^^), which are given by equations ([55]) and (|88p . In 
the equations for the p's and Xp's we should substitute 
jBlp for gi, and 7 for 2e/h/S.^, respectively. 

The macroscopic magnetic moment of a sample with 
N spin-1/2 particles is as follows: 



According to (|83|) the low frequency response of Mz{t) 



M 



(LF), 



it) = ^ [t] cos 29 + p^l^ it) sin 26l) (104) 



where p^'^\t) and p^l' [t) can be expressed in 

terms of the real (%') and imaginary (x") 

parts of their corresponding susceptibilities: 

p(i)(t) = Blf {x'p{^)cofiUjt-x"p{uj)s\nut), and 

^ cos _ (^-j gij^ ^ 

In the same way we obtain from (|88p the low frequency 
response for ^^^^ : 



dM 



{LF) 



dt 



z— !-r!flSin20pL''(t) - -^7^ 

Z 1 1 



(105) 



where p''_^\t) = iB^p ^Xp_ {'^) sinwi + Xp_ (t^) cosa;<j 

The corresponding susceptibilities are given by Eqs. 
(pTjl . ([^ . and or for small high frequency detuning 
((5 < rjfl) by Eqs.dlll), dM]), and ^ where we should 
substitute 7 for 2e/h/S.e- 

The current state of the art allows one to detect low 
frequency oscillations either of Mz with the help of super- 
conducting quantum interferernce devices (SQUIDs)^i2i 
or of dMz/dt by a high quality resonant tank circuit^^. 

(LF) 

Below we write down the explicit form of {t) for 

the case of small high frequency detuning (6 <C ^r)- In 
this case the susceptibilities are as follows: 



1^1 



(0) 



Xp+(w) = -(57 



(0) 

P+ 
T^d{uo) 



1 



(106) 



(107) 



Mz = ^{(Tz) 



(102) 



( ST) 

The stationary magnetization is given by the 



stationary solution of ([88|) or by the stationary part 
of (|89p. By making of the substitutions = I/T2, 



expression^: 



f ST) 

hjBo, we obtain for a well known 



Xp- i^) = i^l 



where d{uj) is defined in (pS]). flu — ■yBi, and 



(108) 



p(0) = Ef_tanh ( (109) 



1 + iT2sy 



where 



(103) 



For Mf^\t) we obtain: 



ALF),,s _ Nh'y 



Mr'{t)^'-^{jBLp)S:^p^^^x 



p2 



[Az(w)coswt-|-B2(w)sinwi] (110) 
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where 



Az{io) = 







-2r2 












2 






2 



(111) 

















' + (r^ + 


(112) 



— ri^ + Tip/Ti. 

VIII. CONCLUSION 

In this paper in the frame of the dressed state ap- 
proach we have analyzed the interaction of a dissipative 
two level quantum system with high and low frequency 
excitations. We have found a linear response of the cou- 

I 



pled TLS and associated photon field system to a signal 
whose frequency is of the order of the Rabi frequency. 
The response of the system exhibits an undamped low 
frequency oscillation, whose amplitude has a clear reso- 
nance at the Rabi frequency with the width being depen- 
dent on the damping rates of the system. The explicit 
expressions for low frequency susceptibility of the TLS, 
which we obtained in the paper, have a rather general 
nature and can be applied for the investigation of spe- 
cific physical systems. As an illustration of our approach 
we study in detail two rather different physical systems: 
the flux qubit and real spin under NMR conditions. 
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Appendix 



A. Calculation of matrix elements of L in the dressed state basis 

With the aid of ®, ^ and (HID, we obtain for {I, N\L\l, N): 



(l,7V|L|l,iV) = 

sir? e{a,N\L\a,N) + cos^ e{h,N - l\L\h,N - I) + sin 6* cos 6* {a,N\L\h,N - 1) + {h,N - l\L\a,N) = 

sm^e[~V^{a,N\a\a,N) +Ti{h,N\a\h, N)] + cos^ 9 [-T i{h,N - l\a\h,N ~ 1) +T^{a,N - l\(j\a,N - 1)] - 

r^sin6lcos6l[(a,7V|cr|6,iV- 1) + {b, N - l\a\a, N)] (113) 



Further transformation requires the substitution of the uncoupled states in (|113p with the dressed states by using 
Eqs. (fTT|) . and p2|) . As a result we obtain: 



{l,N\L\l,N) = 

- (l,iV|cr|l, AT) [Fi sin"* 61 -I- Fi cos'^ 6* -f- sin^ 61 cos^ 6*] - (2, iV|cr|2, TV) sin^ 6* cos^ 61 [F| + F| - 2F^] + 
[(1, 7V|ct|2, N) + (2, 7V|cr|l, N)] sin 9 cos 9 [F| sin^ 6* - F^ cos^ 61 + F<^ cos 26*] + sin^ 9 cos^ 9T ^{l, N + N + l) + 

sni^9cos^ 9T^{l,N - 1|ct|1,7V- 1) + sin"* 6IF4 (2, iV + l|cr|2, TV + 1) +008"* 6IF|(2, A^- 1|ct|2,7V- 1) + 
[(1, N + l|cr|2, TV + 1) + (2, iV + l|cr|l, iV + 1)] Fx sin^ 6'cos6'- 
[(l, iV - l|cr|2, iV - 1) + (2,iV- l|a|l,iV- 1)]F| sin 6* 008^61 (114) 



We do not explicitly write here the other matrix elements of L, such as (2, A^|L|2, TV), (1, A^|L|2, TV), (1, iV|i|l, - 1), 
(2, Af|L|2,iV- 1), (l,iV|L|2,iV- 1), (2, A^|L|1, - 1), which can be obtained by similar procedure. 
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B. Calculations of the averages of Pauli spin operators {(tx), {<^y), and {(Jz) in the dressed state basis 



{az) = Tr {ggz) - ^(1, N\ggz\\, N) + ^(2, N\ggz\2, N) = ^(1, iV|c7|l, N) (1, N\az\l, N) + 

N N N 

^(l,7V|a|2,iV)(2,7V|az|l,iV) +^(2,iV|a|l,7V)(l,iV|az|2,Ar) +^(2,Ar|^|2,iV)(2,7V|a;j|2,7V) (115) 



N N N 



(a,) = Tr {gg,) = ^(1, N\Ga^\l,N) + ^(2, N\aG,\2, N) = 

N N 

^(1, 7V|a|l, iV + 1) (1, TV + l|a,|l, iV> + ^(1, iV|a|l, TV - 1)(1, TV - l|a,|l, 7V) + 

N N 

^(1, TV|a|2, iV + 1) (2, TV + l|a,|l, TV) + ^(1, iV|a|2, TV - 1)(2, iV - l|a,|l, iV) + 

N N 

^(2, TV|a|l, TV + 1) (1, TV + l|a,|2, TV) + ^(2, iV|a|l, TV - 1)(1, TV - l|a,|2, TV) + 

JV JV 

^(2,TV|a|2,TV+l>(2,TV+l|a,|2,TV)+^(2,TV|a|2,TV + l>(2,A^+l|c7,|2,TV) (116) 



N N 

where i — X,Y. 
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